Constraints on Dark Energy Parameters from Correlations of CMB with LSS 



O 

(N 

< 



Hong Li"'^ and Jun-Qing Xia'^ 
^Institute of High Energy Physics, Chinese Academy of Science, P.O.Box 918-4, Beijing 100049, P.R.China 
''Theoretical Physics Center for Science Facilities (TPCSF), Chinese Academy of Science, P.R.China and 
Scuola Internazionale Superiore di Studi Avanzati, Via Beirut 2-4, 1-34014 Trieste, Italy 

In this paper, we combine the the latest observational data, including the WMAP five-year data 
(WMAP5), the baryon acoustic oscillations (BAO) and type la supernovae (SN) "union" compila- 
tion, and use the Markov Chain Monte Carlo method to determine the dark energy parameters. We 
pay particular attention to the Integrated Sache- Wolfe (ISW) data from the cross-correlations of 
cosmic microwave background (CMB) and large scale structure (LSS). In the ACDM model, we find 
that the ISW data, as a complement to the WMAP data, could significantly improve the constraint 
of curvature Qk- We also check the improvement of constraints from the new prior on the Hubble 
constant and find this new prior could improve the constraint of Qk by a factor of 2. Finally, we 
study the dynamical evolving EoS of dark energy from the current observational data. Based on the 
dynamical dark energy model, parameterizing as 'w{a) = wo -\- Wa{l — a), we find that the ACDM 
model remains a good fit to the current data. When taking into account the ISW data, the error 
bars of wo and Wa could be shrunk slightly. Current constraints on the dynamical dark energy 
model are not conclusive. The future precision measurements are needed. 
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I. INTRODUCTION 

Unveiling the origin of the current accelerating expan- 
sion of our Universe is a big challenge for the modern cos- 
mology either theoretically or observationally. The origin 
source which drives the expansion could be attributed to 
a mysterious budget, dark energy. Thus, the nature of 
dark energy is one of the biggest unsolved problems in 
modern physics and has been extensively investigated in 
recent years. 

The measurements of CMB [1-3], LSS surveys [1, [i| 
and SN 0| have provided a lot of high quality data 
at present. These data have been widely used to con- 
strain various cosmological models. However, one should 
keep in mind that the degeneracies of cosmological pa- 
rameters generally exist in almost all cosmological obser- 
vations, i.e., they are not sensitive to single parameters 
but to some specific combinations of them. These de- 
generacies could weaken constraints on the cosmological 
parameters. It is therefore highly necessary to combine 
different probes to break parameter degeneracies so as to 
achieve tight constraints. Furthermore, different obser- 
vations are affected by different systematic errors, and it 
is thus helpful to reduce potential biases by combining 
different probes. 

One of the useful complementary probe is the late- 
time ISW effect [1]. This ISW effect is produced by the 
CMB photons passing through the time-evolving grav- 
itational potential well, when dark energy or curvature 
becomes important at later times. Therefore, the ISW 
effect provides a promising probe for studying the accel- 
eration mechanism of our universe, especially for the dark 
energy and the curvature of Universe. Cross correlating 
CMB with tracers of LSS surveys for detecting the ISW 
effe ct t9|] has been widely investigated in the literature 

In this paper, we will present the constraints on var- 
ious cosmological models from the current observations. 



including the WMAP5 data, SN "union" compilation, 
and recently released BAO data from SDSS DR7, as well 
as the ISW data. The structured of the paper is as fol- 
lowing: in Section II we describe the method and the 
data sets; in Section HI we present our numerical results 
and discussion; finally we give a summary and outlook in 
Section VI. 



II. METHOD AND DATA 

For the parametrizations of dark energy models, we 
adopt ACDM model and dynamical dark energy model 
with CPL parametrization [27| as following: 



WDE(a) = Wo + Wa{l - a) 



(1) 



where a = l/(l + z) is the scale factor and Wa character- 
izes the "running" of the equation-of-state (EoS) of dark 
energy. In the ACDM model, wo — —1 and Wa = 0. 

For dark energy models whose EoS is not equal to 
— 1 during the evolution of Universe, the perturbation 
of dark energy will inevitably exist. The perturbation of 
dark energy has no effect on the geometric constraints. 
However, when including the CMB and LSS data, the 
perturbation of dark energy should be fully considered, 
because the late time ISW effect will be different signif- 
icantly, and will take an important signature on large 
angular scales of CMB and the matter power spectrum 
of LSS . For quintessence- like or phantom-like mod- 
els, in which w does not cross the cosmological constant 
boundary, the perturbation of dark energy is well defined. 
However, when w crosses —1, one is encountered with the 
divergence problem for perturbations of dark energy at 
w = —I. In order to solve this problem in the global 
analysis, we introduce a small positive constant e to di- 
vide the full range of the allowed value of the EoS w into 
three parts: l)?i;>— l-|-e;2) —l — e<w<—l + e; and 
3) u; < -1-e. 
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FIG. 1: One dimensional distributions of Qk and JIa from different data combinations: WMAP5 (balck solid lines), 
WMAP5+ISW (red dashed lines), WMAP5+HST (blue dash-dotted lines), and WMAP5+HST+ISW (purple dotted lines). 



Working in the conformal Newtonian gauge, the per- 
turbations of DE can be described by 

S = -{l + w){e -S^) -m{cl ~w)S , (2) 

II) S 

e - -n{i-3w)e-- 9 + k\-^ + ^) .(3) 

Neglecting the entropy perturbation, for the regions 1) 
and 3), the EOS does not across —1 and the perturba- 
tion is weh defined by solving Eqs. (|2|3p . For the case 2), 
the perturbation of energy density 5 and divergence of 
velocity, 9, and the derivatives of 5 and 9 are finite and 
continuous for the realistic dark energy models. However 
for the perturbations of the parameterizations, there is 
clearly a divergence. In our analysis for such a regime, 
we match the perturbations in region 2) to the regions 
1) and 3) at the boundary and set (5 = and ^ = 0. 
In our numerical calculations we limit the range to be 
|Aw = e| < 10-4 and find our method to be a very 
good approximation to the multi-field dark energy model. 
More detailed treatments can be found in Ref . [28[ . 

The publicly available Markov Chain Monte Carlo 
(MCMC) package CosmoMC^fll] is employed in our 
global fitting, and modifications have been made to in- 
clude dark energy perturbations, and to suit the dark 
energy models which we study [28l [30| . Furthermore, we 
assume purely adiabatic initial conditions in our calcula- 
tions. 

Our most general theory parameter space vector is: 

P = {uJb,t^c,'ds,T,Wl3,Wa,^k,ns,As,cl), (4) 

where uji, = VLi,h? and = ^ch^, in which fif, and ftc 
are the physical baryon and cold dark matter densities 
relative to the critical density, Qk is the spatial curva- 
ture and satisfies flk + flm + ^^A ~ 1, Os is the ratio 



(multiplied by 100) of the sound horizon to the angular 
diameter distance at decoupling, r is the optical depth 
to re-ionization, wq and Wa are the parameters of dark 
energy EoS given by Eq.([T]), and characterize the 
power spectrum of primordial scalar perturbations, is 
the sound speed of dark energy. For the pivot scale we 
set kso = O.OSMpc"^ 

In the computation of CMB we have included the 
WMAP5 temperature and polarization power spectra 
with the routine for computing the likelihood supplied 
by the WMAP team^. We also combine the distance 
measurements from BAO and SNIa. For the BAO infor- 
mation, we use the newly released gaussian priors on the 
distance ratios, rjDyiz) = 0.1905 ± 0.0061 at z = 0.2 
and rs/Dy{z) = 0.1097 ± 0.0036 at z = 0.35, with a cor- 
relation coefficient of 0.337, which were measured from 
the power spectrum for the distribution of the spectro- 
scopic Sloan Digital Sky Survey (SDSS) Data Release 7 
(DR7) galaxy sample and the 2-degree Field Galaxy Red- 
shift Survey (2dFGRS) data[3i]. For the calculation of 
the likelihood from supernova, we use are the "Union" 
compilation (307 sample) j0] and have marginalized over 
the nuisance parameter [31,]. 

For the ISW data, we have included the package 
for calculating the ISW likelihood function provided by 
Ref. [13, which contain a 3.7 cr detection of ISW by cross 
correlating WMAP data with the LSS data sets of 2- 
Micron All Sky Survey sample, SDSS photometric Lu- 
minous Red Galaxies, SDSS photometric quasars, and 
NRAO VLA Sky Survey radio sources. In our calcu- 
lations, we have taken the total likelihood to be the 
products of the separate likelihoods C of CMB, BAO, 
SN and ISW. Furthermore, we make use of the newly 
released HST prior on the Hubble constant, which is 



^ Available at; |http:/ /cosmologist.info/cosmomc^ 



^ Available at the LAMBDA website: 

|http://lambda.gsfc. nasa.gov/ 1 



3 



the measurement of the Hubble parameter by the Near 
Infrared Camera and Multi-Object Spectrometer (NIC- 
MOS) Camera 2 of the Hubble Space Telescope (HST) 
and give Hq = lOOh km s~^ Mpc~^ by a Gaussian like- 
lihood function centered around Hq = 74.2 and with a 
standard deviation a = 3.8 !32|. 



III. NUMERICAL RESULTS 

In this section we present our global fitting results of 
the cosmological parameters determined from the latest 
observational data. 

In Table I we list the constraints on the dark en- 
ergy density ft^ and the curvature ilk in the ACDM 
model from different data combinations. Due to the 
well-known degeneracy between Ha and flk, we obtain 
a weak constraint on the curvature ftk from WMAP5 
data only, flk = -0.036 ± 0.043 (68%C.L.), as shown 
in Fig[TJ Our universe is very close to flatness, which 
is consistent with the prediction of inflation paradigm. 
However, this degeneracy could be broken by adding 
other different cosmological data, such as the large scale 
structure and supernovae data. When we add the ISW 
data into the analysis, we can find that the combined 
constraint from WMAP5-I-ISW is significantly improved 
over using WMAP5 alone, namely the 68% interval is 
rjfc = -0.002 ± 0.016. The constraints of flk and r^A 
improve by a factor of 2.6 and 2.2, respectively. The 
ISW data give the remarkable complementary effect on 
WMAP5 data in constraining cosmological parameters 
and break the degeneracy between ilfc and JIa. 

In Fig H] we also plot the one dimensional posterior dis- 
tribution of flk and Qa from WMAP5-FHST prior. Here, 
we consider the new released HST prior, h = 0.742±0.038 
(1 cr). We find that this new HST prior could play an im- 
portant role in constraining cosmological parameter and 
give much tighter constraints on cosmological parame- 
ters, flk = 0.002 ± 0.007 (Icr). The error bar has been 
shrunk by a factor of 6.6 when compared to the constraint 
from WMAP5 alone. The 95% intervals of fifc and J^a are 
-0.012 < flk < 0.014 and 0.70 < flA < 0.80. We also do 
the calculations using the old HST prior, h = 0.72 ± 0.08 
[33|. Using WMAP5+HST old prior, we obtain that the 
2 a constraints on flk and Ha are —0.045 < ilk < 0.015 
and 0.59 < Ha < 0.80, respectively. The constraints sig- 
nificantly improve by a factor of 2. Furthermore, when 
combining WMAP5, ISW and HST together, the con- 
straint on ilk will improve further: il^ — 0.003 ± 0.006 
at 68% confidence level. 

In Figl2]we present the two dimensional contours of Hk 
and r^A from different data combinations in the ACDM 
model. We find that ftk and Ha are highly correlated 
with each other. The ISW data and the new HST prior 
are very helpful in breaking such degeneracy. 

We now present the constraints on dark energy pa- 
rameters. In Ref . [23| . they considered the constraints on 
dark energy model with constant EoS, and found that 




FIG. 2: Two dimensional contours of Qk and JIa 
from different data combinations within the ACDM 
model: WMAP5 (balck solid lines), WMAP5-HSW (red 
dashed lines), WMAP5-I-HST (blue dash-dotted lines), and 
WMAP5+HST+ISW (purple dotted lines). 
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FIG. 3: Two dimension contours of the EoS parameters of 
dark energy wq and Wa- The black solid lines are given by 
WMAP5+SN+BAO and the red dashed lines are given by 
WMAP5+SN-HBAO+ISW, respectively. 



the ISW data could modestly improve constraint of dark 
energy EoS. As we know, the variation of dark energy 
EoS could also affect the evolution of gravitational po- 
tential, and the ISW effect consequently. It would be 
worth checking the capabilities of ISW data on the con- 
straint of time- varying dark energy EoS. In FigOwe plot 
the two dimensional contours of wq and Wa from differ- 
ent data combinations. The black solid lines are given 
by fitting with the combination of WMAP5 temperature 
and polarization power spectra, supernovae data, as well 
as the newly released BAO data. Due to the limits of the 
precisions of observational data, the variance of wq and 
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TABLE I. Mean 1 a constraints on fifc and Q,a from different data combinations in the ACDM model. 





WMAP only 


WMAP+ISW 


WMAP+HST 


WMAP+ISW+HST 


100 X Q.k 


-3.57 ±4.28 


-0.182 ± 1.61 


0.226 ± 0.65 


0.304 ± 0.630 




0.635 ±0.130 


0.742 ± 0.0566 


0.752 ±0.0243 


0.756 ±0.0228 



TABLE IL Mean 1 a constraints on wq, Wa and f^A from 
different data combinations in the dynamical dark energy 
model. 





WMAP±SN±BAO 


WMAP±SN±BAO±ISW 


Wo 


-1.002 ± 0.136 


-0.970 ±0.130 


■Wa 


0.0096 ± 0.64 


-0.134 ±0.612 




0.732 ± 0.0149 


0.734 ± 0.0144 



Wa are still large, namely the 95% constraints on wq and 
Wa are -1.25 < Wo < -0.72 and -1.47 < Wa < 1.00 re- 
spectively, which are listed in Table II. This result implies 
that the dynamical dark energy models are not excluded 
and the current data cannot distinguish different dark 
energy models decisively. The ACDM model, however, is 
still a good fit right now. 

We also include the ISW data into the calculations. 
From the red dashed lines in Fig|31 we can find that the 
constraints on the dark energy EoS parameters slightly 
improve from WMAP5+SN+BA0+ISW combination. 
The 95% intervals are -1.21 < wq < -0.71 and -1.52 < 
Wa < 0.86. The ISW data help tighten the lower limit 
on Wq and the upper limit on Wa, but not significantly. 
These is due to the constraining power of SNIa and 
BAO, while the current ISW data are still not constrain- 
ing enough. However, the ISW data could be a useful 
probe for constraining other dark energy models, such as 
the early dark energy model [1^. And they also found 
that the future ISW data can improve the constraints on 
the most important cosmological parameters by a factor 

- 1.5. 

Again, we also compare the constraints on dark energy 
parameters between from the new HST prior and from 
the old one. When we combine WMAP5, SN, BAO, as 
well as the old HST prior, we obtain that the 2 a con- 
straints on Wq and Wa are —1.23 < wq < —0.7 and 

— 1.49 < Wa < 0.94, which is slightly weaker than those 
from the new HST prior. 

If dark energy is not the cosmological constant, we 
should consider the dark energy perturbations in our 
calculations. In the framework of the linear perturba- 
tions theory, besides the EoS of dark energy, the dark 
energy perturbations can also be characterized by the 
sound speed, = Spa/ S pa- The sound speed of dark 
energy affects the evolution of perturbations, and leaves 
the signatures on the CMB power spectrum (ss], [H, [11] . 
Thus, in the literature constraining on the sound speed 
cj. from different observational data has been widely in- 
vestigated (e.g. Refs.di,!!^,!!!!). 
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FIG. 4: One dimensional constraints on the sound speed 
from WMAP5±SN±BAO±ISW combination in the dynami- 
cal dark energy model. 



Here, we also constrain the dark energy sound speed 
from current data in the dynamical dark energy model. 
In FigUl we show the one dimensional constraint on 
the dark energy sound speed from the combination 
of WMAP5 temperature and polarization power spectra, 
SNIa and BAO, as well as the ISW data. We can find 
that the constraint on the dark energy sound speed 
are still very weak. The current observational data are 
still not accurate enough. 



IV. SUMMARY 

In this paper, we study the constraints on cosmolog- 
ical parameters from the recently released CMB, BAO 
and SNIa data. Here, we pay particular attention to 
the current ISW data which is the cross correlations of 
CMB with LSS surveys. In the ACDM model, the ISW 
data and the new HST prior are very helpful to break 
the degeneracy between flk and ^Ia- The constraints on 
ilk and significantly improve, when compared to the 
constraints from WMAP 5 alone. 

More importantly, we consider the constraints on the 
dark energy parameters in the CPL dark energy model. 
Here, we fully include the perturbations of dark energy. 
This result implies that the dynamical dark energy mod- 
els are not excluded, the ACDM model, however, is still 
a good fit right now. We find that the ISW data could 
give slight improvement of the constraints on CPL dark 
energy model. But this does not mean that the ISW data 
can not constrain the dynamical dark energy models effi- 
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ciently. Actually, it should be useful for constraining the 
dynamical dark energy models whose EoS deviate 
from the cosmological constant boundary obviously, for 
example, the early dark energy models [25'|, or the dark 
energy model with its EoS w{z) transits sharply during 
its evolution. Besides the dark energy models, the ISW 
data could be also helpful for testing the modified gravity 
theories 40-4^, massive neutrinos [44], the primordial 
non-gaussianity [H, , and so on. 

Furthermore, we compare the constraints between 
from the new HST prior and from the old one. One can 
see that the new HST prior gives the tighter constraints 
on the cosmological parameters. Finally, we check the 
capability of current observational data to constrain the 
dark energy sound speed c^. We find that the sound 
speed is weakly constrained by current observations, and 
thus futuristic precision measurements of the CMB on a 



very large angular scale (low multipoles) are necessary. 



Acknowledgements 

We acknowledge the use of the Legacy Archive for Mi- 
crowave Background Data Analysis (LAMBDA). Sup- 
port for LAMBDA is provided by the NASA Office of 
Space Science. Our numerical analysis was performed 
on the MagicCube of Shanghai Supercomputer Center 
(SSC). We thank Yi-Fu Cai for helpful discussions. This 
work is supported in part by the National Natural Science 
Foundation of China under Grant No. 10803001, 973 
Program under Grant No.2010CB833000 and the Youth 
Foundation of the institute of high energy physics under 
Grant Nos. H95461N. 



E. Komatsu et al. [WMAP Collaboration], Astrophys. J. 
Suppl. 180, 330 (2009). 

J. Dunkley et al. [WMAP Collaboration], Astrophys. J. 
Suppl. 180, 306 (2009). 

G. Hinshaw et al. Astrophys. J. Suppl. 180, 225 (2009). 
M. Tegmark, et al, Astrophys. J. 606, 702 (2004); 
M. Tegmark, et al, Phys. Rev. D 74, 123507 (2006). 
S. Cole, et al, Mon. Not. R. Astron. Soc. 362, 505 (2005). 
M. Kowalski et al Astrophys. J. 686, 749 (2008). 
M. Hicken et al, Astrophys. J. 700, 1097 (2009). 
R. K. Sachs and A. M. Wolfe, Astrophys. J. 147, 73 
(1967). 

R. G. Crittenden and N. Turok, Phys. Rev. Lett. 76, 575 
(1996) 'arXiv:astro-ph/9510072\ 

S. P. Boughn, R. G. Crittenden and N. G. Turok, New 
Astron. 3, 275 (1998) arXiv:astro-ph/9704043 . 
M. Zal darriaga and U. S eljak, Phys. Rev. D 59, 123507 
(1999) [arXiv:astro-ph79 810257 . 

W. Hu and T. Okamoto, Astrophys. J. 574, 566 (2002) 
farXiv:astro-ph/0111606 . 

Y. S. Song, A. Cooray, L. Knox and M. Zaldarriaga, As- 
trophys. J. 590, 664 (2003) ar Xiv:astro-ph /0209001 . 
W. Hu, Phys. Rev. T>^6W, 02300J (2002) 



arXiv:astro-ph/0108090 



557, L79 (2001) 



W. Hu, Astrophys. 
'arXiv:astro-ph/0105424 . 
N. Afshordi, Y. S. Loh and M. A. Strauss, Phys. Rev. D 
69, 083524 (2004) arXiv:astro-ph/0308260; . 
E. Gaztanaga, M. Manera and T. Multamaki, 
Mon. Not. Roy Astron. Soc. 365, 171 (2006) 
arXiv:astro-ph/0407 022 . 

Vielva, E. Martinez-Gonzalez and M. Tucci, 



'arXiv:astro-ph/0408252 

D. Pietrobon, A. Balbi and D. Marinucci, Phys. Rev. D 
74, 043524 (2006) arXiv:astro-pli/0606£75j. 
J. D. McEwen, P. Vielva, M. P. Hobson, E. Martinez- 
Gonzalez and A. N. Lasenby, Mon. Not. Roy. Astron. 
Soc. 373, 1211 (2007) arXiv:astro-ph/0602398'. 
T. Giannantonio et al, Phys. Rev. D 74, 063520 (2006) 
arXiv:astro-ph/0607572 . 

A. Rassat, K. Land, O. Lahav and F. B. Ab- 



[23; 

[24 
[25; 

[26; 

[27 

[28; 

[29 
[30 

[31 

[32; 

[33 
[34' 
[35; 

[36 

[37 

[38; 
[39; 

[40 



dall a, M;on. No t. Roy Astron. Soc. 377, 1085 (2007) 
arXiv:astro-ph /06 109 11. 

S. Ho, C. Hirata, N. Padmanabhan, U. Seljak 
and N. Bahcall, Phys. Rev. D 78, 043519 (2008) 
tarXiv:0801.0642 [astro-ph]]. 

B. R. Granett, M. C. Neyrinck and L Szapudi, 
arXiv:0805.3695 [astro-ph]. 

J. Q. Xia, M. Viel, C. Baccigalupi and S. Matarrese, 
arXiv:0907.4753 [astro-ph. CO]. 

L. Pogosian, P. S. Corasaniti, C. Stephan-Otto, R. Crit- 
tenden and R. Nichol, Phys. Rev. D 72, 103519 (2005) 
arXiv:astro-ph/0506396 . 

M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10, 
213 (2001); E. V. Linder, Phys. Rev. Lett. 90, 091301 

(2003) . 

G. B. Zhao, J. Q. Xia, M. Li, B. Feng and X. Zhang, Phys. 
Rev. D 72, 123515 (2005) arXiv:astro-ph/0507482 . 
A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002). 
J. Q. Xia, G. B. Zhao, B. Feng, H. Li and X. Zhang, Phys. 
Rev. D 73, 063521 (2006) arXiv:astro-ph/0511625 . 
E. Di Pietro and J. F. Claeskens, Mon. Not. Roy. Astron. 
Soc. 341, 1299 (2003). 

A. G. Riess et al, Astrophys. J. 699, 539 (2009) 
[arXiv:0 905.0695 [astro-ph.CO]]^ 

W. J. Percival et al, arXiv:090 7l660l [astro-ph.CO]. 
W. L. Freedman, et al, Astrophys. J. 553, 47 (2001). 
J. Weller and A. M. Lewis, Mon. Not. Roy. Astron. Soc. 
346, 987 (2003) arXiv;astro-ph/0307104; . 
S. DeDeo, R. R. Caldwell and P. J. Steinhardt, Phys. 
Rev. D 67, 103509 (2003) [Erratum-ibid. D 69, 129902 

(2004) ] arXiv:as tro-ph/030 1284 . 

J. K. Erickson, R. R. Caldwell, P. J. Steinhardt, 

C. Armendariz-Picon and V. F. Mukhanov, Phys. Rev. 
Lett. 88, 121301 (2002) arXiv;astro-ph/0112438 . 

W. Hu and R. Scranton, PhysTRev. D 70, 123002 (2004) 
^arXiv:astro-ph /0408456;^ . 

J. Q. Xia, Y. F. Cai,' T. T. Qiu, G. B. Zhao and 
X. Zhang, Int. J. Mod. Phys. D 17, 1229 (2008) 
[arXiv;astro-ph/0703202 . 

B. Jain and P. Zhang, Phys. Rev. D 78, 063503 (2008) 
,arXiv:0709.2375 [astro-ph]]. 



6 



[41] N. Afshordi, Phys. Rev. D 70, 083536 (2004) [44] J. Lesgourgues, W. Valkenburg and E. Gaztanaga, Phys. 
[arXiv:astro-ph/0401166 . Rev. D 77, 063505 (2008) arXiv:0710.5525 [astro-ph]]. 

[42] 
[43] 



Zhang, Phys. Rev. D 73, 123504 (2006) [45] A. Slosar et ai, JCAP 0808, 031 (2008). 



arXiv:astro-ph/0511218;. [46] J. Q. Xia et ai, arXiv:1003.345T1 

Schmidt, M. Liguori and S. Dodelson, Phys. Rev. D 
76, 083518 (2007) arXiv:0706.1775l [astro-ph]]. 



